Gamma-Ray Burst Afterglow: Polarization and Analytic Light Curves 
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ABSTRACT 

GRB afterglow polarization is discussed. We find an observable, up to ~ 10%, 
polarization, if the magnetic field coherence length grows at about the speed of light 
after the field is generated at the shock front. Detection of a polarized afterglow 
would show that collisionless ultrarelativistic shocks can generate strong large scale 
magnetic fields and confirm the synchrotron afterglow model. Non-detection, at ~ 1% 
level, would imply that either the synchrotron emission model is incorrect, or that 
strong magnetic fields, after they are generated in the shock, somehow manage to stay 
un-dissipated at "microscopic", skin depth, scales. Analytic lightcurves of synchrotron 
emission from an ultrarelativistic self-similar blast wave are obtained for an arbitrary 
electron distribution function, taking into account the effects of synchrotron cooling. 
The peak synchrotron fiux and the fiux at frequencies much smaller than the peak 
frequency are insensitive to the details of the electron distribution function; hence 
their observational determination would provide strong constraints on blast wave 
parameters. 



Subject headings: gamma rays: bursts — magnetic fields — shocks 



1. Introduction 

X-ray, optical and radio emission following gamma-ray bursts (GRBs) are in broad agreement 
with models based on relativistic blast waves at cosmological distances ([Waxman 1997a , Wijers 



Rees fc Meszaros 1997 , Vietri 1997b , Reichart 1997 , Katz &: Piran 1997 , Sari, Piran, &: Narayan 
1998). In these models, the energy released by an explosion, ~ lO^^erg, is converted to kinetic 



energy of a thin baryon shell expanding at ultra-relativistic speed. After producing the GRB, the 
shell impacts on surrounding gas, driving an ultra-relativistic shock into the ambient medium. 
In what follows, we refer to the surrounding gas as interstellar medium (ISM) gas, although 
the gas need not necessarily be inter-stellar. The expanding shock continuously heats fresh gas 
and accelerates relativistic electrons, which produce the observed afterglow radiation through 



synchrotron emission ( Paczyhski &: Rhoads 1993 , Meszaros &: Rees 1997 , Vietri 1997a ). 



To match the observations, the magnetic field behind the shock has to be ~ 10% of 
equipartition with the shock-heated, compressed ISM. What is the origin of this field? The 
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shock-compressed ISM field is many orders of magnitude smaller than needed. The magnetic field 
frozen into the initial GRB fireball loses strength by the time the afterglow stage begins, and it is 
in a wrong place anyway. During the afterglow, the decompressed GRB field is located far behind 
the shock, while most of the energy is in the recently shocked ISM. Therefore, the magnetic field 
most likely must be generated in and by the blast wave. If the coherence length of the generated 
field is comparable to the thickness of the blast wave, the radiation will be polarized. An ~ 10% 
degree of polarization is expected. This is significantly smaller than the maximal synchrotron 
polarization, ~ 70%, because the emitting region is thin and broad; it must be covered by ~ 100 
mutually incoherent patches of magnetic field. 

In a paper on microlensing of GRB afterglows, Loeb & Perna (1998) have mentioned the 
possibility that the afterglows are polarized. Here we estimate the degree of polarization (§4). 
This paper also provides (§3) exact analytic afterglow lightcurves for an arbitrary electron 
distribution function, including the effects of electron cooling. In §2 we describe the underlying 
model assumptions. We discuss the implications of our results to afterglow observations in §5. 
Most of the details of our derivations are given in appendices A-E. 

2. The blast wave model 

A strong blast wave is fully specified by two parameters: the blast wave energy 
E = 10^'^£^52 erg, and the ISM density rij = Ini cm~^. With sufficient accuracy, the unshocked 
ISM may be taken to be cold unmagnetized hydrogen. To calculate the synchrotron emission we 
need to know the fraction of energy in magnetic fields ^s, and in electrons .^ei and the shape of the 
electron distribution function (a function fe{z) with first two moments equal to 1). We include 
Cb, Ce, and feiz) in the list of independent parameters. In principle, these are determined by the 
blast wave energy and the ISM density, but a theory of strong collisionless shocks is not available 
(Sagdeev 1966, Krall 1997). 

The plasma flow in the shocked ISM is assumed to be described by the Blandford-McKee 
( Blandford fc McKee 197^ ) self-similar solution. The Lorentz factor of the shock wave T, the 
Lorentz factor of the flow 7, the proper energy density e, and the proper number density n for all 
space-time points in the shocked plasma are given in Appendix A. 

We assume that magnetic fields and electrons are described by simple scalings. The magnetic 
field is B'^/8tt = £,3^, and S,b is the same in all space-time points. The electron distribution 
function in the local rest frame has the same shape in all space-time points on the shock front, 
after the shock passage it evolves by adiabatic and synchrotron cooling. At the shock front, the 
mean energy of an electron in the local rest frame is 7emeC^ = Cec/n, and is constant. The 
shape of the electron distribution function is not specified at this point. We include feiz) into the 
definition of the synchrotron emission function F. The synchrotron power per unit frequency per 
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electron emitted in the local rest frame is 

V3 e^B ^ , oj 

o 2^ ^~ 



= ir—n-l, (1) 



c..= 52M. (2) 

2meC ' ^ ^ 

but F is not the standard dimensionless synchrotron emission function given in Rybicki & 
Lightman (1979). F depends on the shape of the electron distribution function (B7). 

Given the set of blast wave parameters, we will measure time, frequency, and spectral 
luminosity (energy per time per frequency) in units of 

T = c-'( -— %1 = 5.5 X 10^ ( ^] s, (3) 
\S-K nimp& J \ ni J 

- - (M' - 3. X m'^ „;/v.. (.) 
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E,^-^ {^] iT{n.rir-E = 2.2 x 103^ (|^) n^E.^erg. (5) 
2V67r \mpj " VO.l/ 



The formal origin of these units is explained in Appendices A, B. Their physical meaning is 
illustrated by the following order of magnitude statement. At observed time to/T = 1, the blast 
wave slows down to Lorentz factor 2; it radiates at frequency oo/ijOq = 1, with spectral luminosity 
L/Eq = 1. Our analysis is restricted to the ultrarelativistic stage, that is to dimensionless observed 
times to 1- 



3. Lightcurves 

We first calculate in §3.1 synchrotron emission of the blast wave neglecting radiative cooling 
of electrons, i.e. assuming that the shape fe{z) is the same in the entire shocked plasma. In §3.2 
we relax this assumption: fe{z) is determined at the shock front and evolves by synchrotron and 
adiabatic cooling thereafter. 

Our analytic lightcurves are exact under the following assumptions: (i) The blast-wave 
hydrodynamics is described by the Blandford-McKee self-similar solution ( piandford &: McKee 



19761) ; (ii) The magnetic field energy density is a fixed fraction of the total energy density. 



independent of space and time; (iii) The electron distribution function is determined at the shock 



front and evolves afterwards only through adiabatic and synchrotron cooling. Granot, Piran, & 



Sari 1998 numerically derived exact lightcurves for power-law electron distribution, under the 
assumptions described above and neglecting electron cooling. It should be emphasized that since a 
theory of strong collisionless shocks is not available at present, none of the above assumptions can 
be justified. Thus, the numerical values (e.g. of the peak flux and peak frequency as function of 
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time) derived here under these assumptions are not necessarily more accurate than those obtained 
by order of magnitude estimates (e.g. Waxman 1997b| , Sari, Piran, & Narayan 1998, Wijers & 
Galama 199"^ ). 



The exact analytic lightcurves are useful because they allow us to determine which afterglow 
characteristics are strongly dependent on the details of the electron distribution function, and 
which are insensitive to these details and depend mainly on the global blast wave parameters (i.e 
the blast wave energy, the ambient medium density, and the energy fractions carried by electrons 
and magnetic fields). 



3.1. Adiabatic Lightcurve 

In Appendix B, we derive an expression for the spectral luminosity, neglecting synchrotron 
cooling of electrons. At observed time to after the gamma-burst, at frequency iv, distant 
observer (with negligible redshift) infers a selfsimilar narrow-band luminosity of the blast wave 
Luj{to) = LA{iotl^'^), where 

La{uj) = 48 r da o^(l + 7a'^y'^F[ 2a(l + 7a^)uj ]. (6) 
Jo 

We show the lightcurves in Fig.l. The three curves correspond to different doubly normalized 
electron distribution functions: 

1. Power law, index p = 2.4: f^ = fp = 31.2^2 (1 + 122zP+2)-i 

2. Maxwelhan: = /a/ = IS.bz'^e'^' 

3. Mixed: fe = 0.7/m + 0.3/p 

Note, that our "power-law" distribution, for which /g oc for z ^ 1, includes a low energy tail, 
/e ~ z'^ for z <^ 1. We believe the inclusion of such a "thermal" low energy tail is more realistic 
than assuming a sharp cutoff of the electron distribution below a certain minimum z value. 



3.2. Lightcurve with synchrotron cooling 

At early times or at high frequencies, synchrotron cooling of the electron distribution 
function will have a noticeable effect on the lightcurve. We calculate the nonadiabatic lightcurve, 
Lu){to) = L]\[A{i^,to) in Appendix C, neglecting effects of cooling on the blast wave propagation: 

LNA{uj,to) = 192 r dy y^ da a\l + 7a^y^ [ dzo fe{zo) i^o[ 2a(l + 7a^)u:tl'^ z~^ ], (7) 
Jo Jo J 
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where 

z-^ = z^' + A{8to)-'/^a-'y-\l - y''/^), (8) 

and 

Scaled spectra at different observed time to are shown in Fig. 2 for the power-law (and in Fig. 
3 for the mixed) electron distribution function of §3.1 and A = 0.01. At high frequencies, eqs. 
(7), (8) predict a power- law luminosity L oc uj~p/^ for an electron distribution function with a 
power-law tail of index p. 



3.3. Observables 

From Eqs. (3)- (5) and Figs. 2, 3, an afterglow at redshift zi,, observed tjay days after the 
7-burst, will show maximal flux at a frequency 

~ 3 X 10-^ (It) ' (si) £5=''*i''H- 

The maximal flux does not depend on the time of observation. In a fiat universe with Hubble 
constant Hq = 75 km/s/Mps, 

= ^ — ; — TTT "^1 -^52mJy. (11) 



^/2-l J VO.l 

As seen in Figs. 1, 2, 3, the peak flux F^^jn is robust, i.e. it is independent of the details of the 
electron distribution function. The flux below the peak is also robust, and for i/ ^ i/^ it is given 

by 

P r.J V^^b-l V'fl + Zb\-'/'fCeY^^^f^B\'/' 1/2^5/6,1/2 1/3 , 

The peak frequency is model-dependent, and may differ by an order of magnitude between 
different electron distribution functions with similar ^g- The spectral shape (equivalently the 
time profile) above the peak is also strongly dependent on the details of the electron distribution 
function. 



4. Polarization 

Synchrotron radiation is highly polarized (Ginzburg 1989), but for this polarization to be 
measurable in an unresolved source, the magnetic field coherence length should be comparable to 
the source size. Here we show that GRB afterglows might be polarized. An ~ 10% polarization 
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seems to be an upper bound, corresponding to a coherence length that grows at about the speed 
of hght after the field is generated at the shock front. 

Qualitatively, our polarization analysis can be summarized as follows. Suppose that the 

magnetic field coherence length in the local rest frame is I ~ cr, where r is the proper time after 
the shock. The extension of the emitting region transverse to the line of site is ~ 5cr. There are 
~ 50 coherent patches. The degree of polarization is ~ 60%/\/50 ~ 10%. If the coherence length 
is smaller than the proper time, I ^ ecr, e < 1, the degree of polarization is decreased to 

n ~ lOe^/^ %. (13) 

The degree and direction of polarization should depend on time, the polarization coherence time 
should be ~ etg- 

4.1. Magnetic field generated by a relativistic coUisionless shock 

As far as we know, magnetic field generation in coUisionless shocks is not understood. It seems 
possible that, at the shock front, Weibel instability generates near equipartition (^b ~ 0.1) small 
scale (~ skin depth 6, here 6 = c/up, ujp ^ ne^/^pirip ^ ne'^/^e^e) magnetic fields. By magnetic 
moment conservation, electrons are accelerated to near equipartition energies ( relativistic version 
of Sagdeev 1966, Kazimura et. al. 1998). The ultimate fate of the field many skin depth behind 
the shock front is not clear. What happens to the magnetic field coherence length /, and to the 
magnitude 5 at a distance A 2> behind the shock front? Three scenarios seem to make sense: 

1. The generated field dies out after finishing its job of isotropizing the plasma and bringing it 
to a state given by the shock jump conditions. 

2. The magnitude stays at about equipartition, the coherence length stays at about the skin 
depth. 

3. The magnitude stays at about equipartition, the coherence length grows as Z ~ A. 

Scenario 1 is not consistent with the synchrotron emission model for the afterglow, because 
too little synchrotron radiation is produced by a skin-deep shell with strong magnetic fields. 
Scenario 2 means unpolarized radiation. We will evaluate the degree of polarization for scenario 3. 

4.2. Coherent patch 

Assume that two events in the shocked ISM belong to the same coherent patch if the difference 
in their proper times elapsed after the shock passage St, and their spatial separation transverse to 
the line of site Sh are small: 

St < erT, (14) 
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5h < ehCT, (15) 

e-r, e/i < 1, T is the averaged proper time since the shock. By the proper time of an event in the 
shocked ISM we mean the proper time after the shock of a fluid particle at this event. 



4.3. Degree of Polarization 

As shown in Appendix D, for the emission event with observed time to, the proper time r (in 
units of T) and the transverse distance (in units of cT) are 

T=1.15t^/8a5/V/'(l-y9/4), (16) 
/i = 1.83t^/8al/4(l_„2)l/2yl/4_ 

Here the meaning of dimensionless variables a, y is unimportant, what matters is that the 
luminosity is given by the integral (B20) over a and y: 

L{uj) = 192 dy y^ da a^{l + 7a'^)-'^F[ 2a(l + 7a'^)uj ]. (18) 
Jo Jo 

Now we can separate the full luminosity (18) into coherent parts according to the criterion 
(14), (15). This gives, approximately, the degree of polarization. By numerical simulations 
(Appendix E), 

n ~ lOe/^e^^ %. (19) 



5. Summary of results 
5.1. Light curves 

We have derived simple, exact analytic afterglow lightcurves for an arbitrary electron 
distribution function, including the effects of electron synchrotron cooling, Eq. (7), and in the limit 
where synchrotron cooling is negligible, Eq. (6). Our lightcurves are exact under the following 
assumptions: (i) The blast-wave hydrodynamics is described by the Blandford-McKee self-similar 
solution ( piandford &: McKee 197^ ); (ii) The magnetic field energy density is a fixed fraction of 
the total energy density, independent of space and time; (iii) The electron distribution function 
is determined at the shock front and evolves afterwards only through adiabatic and synchrotron 
cooling. 

We have shown, see. Fig. 1, 2, 3, that the peak synchrotron flux, Eq. (11), and the 
synchrotron flux at frequencies well below the peak flux, Eq. (12), are insensitive to the details 
of the electron distribution function. Since the peak flux is also insensitive to the details of 
the blast wave hydrodynamic profiles (Fig. 1 in Waxman 1997c| ), the peak flux and the flux at 
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frequencies well below the peak depend mainly on the global blast wave parameters: blast wave 
energy, ambient medium density, magnetic field and electron energy fractions [cf. Eqs. (11,12)]. 
Observational determination of these fluxes would therefore provide strong constraints on blast 
wave parameters. The numerical value of the peak flux derived here, Eq. (11), is similar to that 



derived in Granot, Piran, & Sari 1998, within a factor ~ 3 of the values given in Waxman 1997b 



and [Wijers fc Galama 1998 , and a factor ~ 10 smaller than given in Sari, Piran, &: Narayan 1998| . 



The discrepancy with pari, Piran, &: Narayan 1998| is mainly due to the fact that it is assumed by 
these authors that the spectral width of the observed spectrum at fixed time is determined mainly 
by the intrinsic spectral width of synchrotron emission, while the actual width is significantly 
larger and dominated by contribution to the observed spectrum at given time from different 
space-time points with different plasma conditions. 

The peak frequency and spectral shape at frequencies above the peak are strongly dependent 
on details of the electron distribution function, see Figs. 1, 2, 3. Furthermore, the peak frequency 
is also strongly dependent on the details of the blast wave hydrodynamic profiles (Fig. 1 in 
Waxman 1997c). This, and the fact that the spectral peak is flat, imply that observational 



determination of the peak frequency and of spectral features above the peak at a given time can 
not be used directly to constrain global blast wave parameters. These features would mostly 
provide information on the electron distribution function. The numeric value of the peak frequency 
derived here, Eq. (10), is within a factor ~ 3 of the values given in Sari, Piran, fc Narayan 1998 



Granot, Piran, fc Sari 1998| and [Wijers &: Galama 1998| , and a factor ~ 10 smaller than given in 
Waxman 1997b . The discrepancy with Waxman 1997b| is due mainly to the fact that it is assumed 



in Waxman 1997b that the spectral peak is close to the synchrotron frequency of electrons with 
average Lorentz factor, while the actual peak is closer to the synchrotron frequency of electrons 
near the peak of the electron distribution function. 

The break frequency (the frequency where the high-frequency spectrum changes the slope 
from (p — l)/2 to p/2 due to synchrotron cooling) is not prominent. The transition to the p/2 
slope occurs in a manner that strongly depends on the details of the electron distribution function. 



5.2. Polarization 

If the observed afterglows are indeed synchrotron emission from ultrarelativistic blast waves 
propagating into ISM, the magnetic field needed to account for the emission must be generated by 
the blast wave. If the coherence length of the generated field grows at about the speed of light 
after the field was generated at the shock front, afterglows should be noticeably polarized. 

We thank John Bahcall for a discussion that initiated this study. Our work was supported by 
NSF PHY-9513835. EW was also supported by the W. M. Keck Foundation. 
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A. Ultrarelativistic blast wave. 

The Blandford & McKee (1976) solution can be deseribed as follows. 

Let (t, r, 0) be the space-time coordinates in the blast frame, 6 is the polar angle which is 
assumed to be small, with ^ = in the observer direction. Let E be the energy of the blast wave, 
and Hi the unshocked ISM number density, c = 1. The shock front propagates into the ISM with 
a Lorentz factor T that decreases with time according to 

.2,3 17 E 



r^r = —33- = (ai) 

Define a similarity variable 



Stt riimp 



X = 8r\l--). (A2) 
The shocked region is % > 1, and the fluid flow in the shocked region is given by 

1' = (A3) 

e = 2r\-^y^^nimp, (A4) 

n = 2V2Tx~^/^ni. (A5) 

Here 7 is the Lorentz factor of the flow, e is the proper energy density, n is the proper number 
density. 

In Appendix B, we use (F, 7,tobs) as independent variables instead of {t,r,9). Here iobs is the 
time at which a photon emitted at (t, r, 9) is observed; with sufficient accuracy, 

tohs = t-r + r— = t-r + t—. (A6) 
The coordinate transformation is (old coordinates in terms of new coordinates) 

t = (A7) 

^ = (i-T^)^ = (i-T^)^r-V3, (A8) 

^obs 1 ^obs p2/3 _ 1 ^^g-j 



2 t 16Y T 167^ 



From (A3) 



X = ^- (AlO) 

With X from (AlO), expressions (A4), (A5) give the dependent quantities, energy density and 
density, in terms of new independent variables. In new coordinates, the space-time domain of the 
shocked fluid (t > 0, x > 1) ^ > 0) is given by 

00 > tobs > 0, (All) 
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'o. \ -3/8 

oo > r > ( ^) , (A12) 



B. Lightcurve 

Here we calculate the lightcurve of synchrotron emission from an ultrarelativistic blast wave. 
The physical assumptions of the model are discussed in the main text. 

It is convenient to start from the following expression for the total emitted energy 

Er = Jr^dr J 2^9(16 j dt n J duj'P{Lo',je,B)^ (Bl) 
The factors in (Bl) are: 

1. Total number of emitting electrons at a given time 

r'^dr27Tede -fn. (B2) 

2. Total energy emitted by one electron 

/ ^ [du;'Piu;'ne,B)^, (B3) 

where to is the photon frequency in the burst frame, and u' is the frequency in the local rest 
frame, 

J = 1^.. (B4) 

P(u;',7e,-B) is the synchrotron radiation spectral power in the local rest frame emitted at 
the frequency u;', by one electron from a distribution with a mean Lorentz factor 7e, in the 
magnetic field B. It is given by (Rybicki & Lightman 1979) 

The dependence of emission on the mean Lorentz factor of electrons is shown explicitly. The 
dependence on the detailed distribution function of electrons is hidden in the definition of 
the synchrotron emission function F{x). Namely, we define 

F{x) = j dzfe{z)Fo{^), (B7) 
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with Fq (x) being the standard synchrotron emission function 

Fo{x) = x i^5/3(0- (B8) 

Jx 

In (B7), the normahzed electron distribution function /g is defined by the following 
expression for the probability for the electron to have a Lorentz factor jei 

dProbability^J_^ /7^\ (B9) 



le 



3. The last factor in (Bl) is the ratio of infinitesimal solid angles in the local rest and blast 
frames: 

dn' _ 

In ~ (1+72^2)2- (BIO) 

We assume that magnetic fields and electrons take up a fixed fraction of the proper energy 
density e{r,t): 

^=^Be, (Bll) 

and 

'yenm.eC^ = (B12) 

We also assume that the normalized electron distribution function fe{z) in the shocked ISM is 
fixed. These assumptions might be approximately correct when synchrotron cooling becomes 
unimportant at later stages of the afterglow. 

We use the Blandford &; McKee (1976) selfsimilar solution for the Lorentz factor 7, density 
n, and energy density e, and change the independent variables in (Bl) from {t,r,9,Li>') to 
(r,7,tobs)'^)- Using Appendix A, we get 

17 E 

Er = — 5 / du; dtohs 

Jo Jo 



48 mnipC^ 



/ / d-f -f-^ nD^P(D-^u;,-fe,B). (B13) 

Here D is the "Doppler factor" 

^=n^=2^(^+2T''^''v)", (B14) 

and T is the characteristic time of the blast wave introduced in Appendix A. 
The spectral lightcurve is defined as luminosity per unit frequency: 

^c.(tobs) = -J^- (B15) 
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Prom (A14), 

Now wc use expression (B5) for the synchrotron power P, (Bll) and (B12) for B and 7e, (A4), 
(A5), (AlO) for e and n. Also, from now on we will denote by to the observed time measured in 
units of T. We also define the frequency and spectral luminosity units, equations (4), (5). These 
are devised to get rid of constant factors in the resulting expression for the luminosity. We denote 
the frequency lh in units of loq by u, and the spectral luminosity L^{tohs) in units of £^0 by L^{to)- 
Then (B16) takes the following form 

Define new integration variables x and y: 7 = {T/\/2)y, T = {^to)~^^^x~^/'^ . We obtain a 
selfsimilar spectral lightcurve 

L^{to) = LAiojtl/^). (B18) 

Here 

• 1 1 2 2 

La{co) = 192 / dx / dy y35/i2(7 + l^y^Fl 2x^y-^''/^\7 + ^)io ]. (B19) 

Jo Jx X X 

With only an ~ 2% error in the resulting luminosity, we can replace the indices 35/12 and 37/12 
by 3, and get a simpler expression 

La{oj) = 192 dy y^ da a^(l + 7a^y^F[ 2a(l + 7a^)u} ], (B20) 
Jo Jo 

where a = x/y. Prom (B20), the adiabatic lightcurve is 

La{oj) = 48 da a^(l + la^)-'^F[ 2a(l + 7a^)uj ]. (B21) 
Jo 

C. Synchrotron cooling 

Synchrotron plus adiabatic cooling of an electron with Lorentz factor ^^l is described by 

c?7ez l-ieldn 4 52 2 in-i\ 

= 3 " 3^^' w:^^^" ^^^^ 

where r is the proper time of the fluid element at the electron's location. We have dr = dt/^. 
Prom (A5), dlnn = dlnP- (5/4)dlnx. Prom (D8), cZlnx = Mint. Then 

d^el I lei 4 2 1 fdlnF 5\-i 



dn 3 n 3 STrmgC^ 777 V 
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Using (Al), (Bll) 

— = -— + —— ^B7^— (C3) 
dn 3 n 39 ^' ^yn 

To integrate, we need to express t, e, and 7 in terms of n. Let 70, uq and cq be Lorentz factor, 
proper density, and proper energy density at the shock passage time to. Prom (Al), (A3)-(A5), 

70 = ^ro, (C4) 

no = 2V2roni, (C5) 

eo = 2Tlnimp(?, (C6) 
where = T^/tg. Prom (D8), (Al), (A3)-(A5), 

7 = lo{t/to)-"\ (C7) 

n = no(tAo)-''/', (C8) 

e = eo(«Ao)-''/'. (C9) 

Now (C3) can be written as 

d^ei l7e; , 8 ar heo ( n \-''"'\2 .p..^ 



dn 3 n 39 mgC 7o?^o V^^o 

and integrated 



- = - + ^^ei.ee^(i-y^n (cii) 

Here z = ^ei/^e, with 7e defined by (B12), y is defined by (Dl). Plug in (C4)-(C6), express Tq in 
terms of a, y, and to- We get 

= z^' + A(8i,)-V2a-iy-2(i _ yi9/6)^ (C12) 

A is defined by (9), a is defined by (D2). With the synchrotron coohng given by (C12), the 
spectral luminosity is 

L^it) = 192 dy y^ da 0^(1 + la^)'^ [ dzofe{zo)Fo[ 2a(l + 7a'^)u;tl/y z'^ ]. (C13) 
Jo Jo J 



D. Transverse distances and proper times 
Adiabatic lightcurve (B20) is an integral over dimensionless variables y and a: 

. = ^. (Dl) 
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ay = (8to)-i/2r-V3_ P2) 

To calculate polarization using (B20), we have to express the distance from the observer - burst 
center line h, and the proper time since the shock passage r, in terms of y and a. This is done 
here. 

The transverse distance is ^ = rO, and from (A8), (A9) 

h = V2rr-2/3(t,r2/3 _ (D3) 

and from (Dl), (D2) 

h=^T{8tof'{ayy/\l-a')'/\ (D4) 



1-^ = 87-7, (D5) 



Now the proper time... Equation of motion of a shocked particle is 

dr 1 r 

It^ 2^^ t 

here t is the burst frame coordinate time. Integration gives r = t — Ct^. Since the shock front is at 

R = t-:^, (D6) 

we get 

where Iq is the burst coordinate time at which the particle was shocked. The similarity variable at 
the particle is 

and 

Using (D8), (AlO), (A7), we get 

r = ^T{8tof/'a'/V/\l - y'/'). (DIO) 



E. Degree of Polarization 

To estimate the degree of polarization 11 we use the adiabatic lightcurve and assume 
F{u>) ~ Lo"^ in (B20). The latter simplification should not lead to a large error, because 11 turns 
out to be approximately the same in the relevant range l>s>— 1/3. With these assumptions, 
the degree of polarization can be estimated as 

n = ^ei£. (El) 
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Here the first factor is polarization of a power-law emission from one patch, L is the total 
unpolarized luminosity, and Cll is the polarized luminosity correlator. Up to an irrelevant factor, 

L = i y daa^-'{l + 7a^)-^-', (E2) 

Cll = J daidyida2dy2{yiy2f{aia2f-'[{l + 7a?)(l + Ta^)]-^-^ C12 min(l, |^). (E3) 

Here the min-term comes from the azimuthal angle integral, r = (ri + T2)/2, h = {hi + /i2)/2. 
C12 is the normalized magnetic field correlator, for which we take a simple form corresponding to 
equations (14), (15). 

C12 = 0{\Ti - T2\ - e^T)e{\hi - h2\ - ehT). (E4) 
We calculated (El) numerically for different values of s, th and e,-. 
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Fig. 1. — Adiabatic lightcurves (6) for different electron distribution functions (§3): power-law P, 
Maxwellian M, and mixed MP. This graph can be interpreted as luminosity at a given frequency 
as a function of time, or as the spectrum at a given time. 
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Fig. 2. — Nonadiabatic lightcurves (7) for A = 0.01, for different observed times, and for a power- 
law electron distribution function (§3). Adiabatic lightcurve is shown for comparison. Nonadiabatic 
curves are marked by the Lorentz factors of the shock front at the time when observed photons 
were emitted from the shock front from ^ = 0. Observed time in days is tday = 80(£'52/ni)^/^r'~^/^. 
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Fig. 3. — Same as Fig. 2, but for the mixed (power-law plus Maxwellian) electron distribution 
function (§3). 



